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^— «, ■ Universal scaling laws form one of the central issues in physics. A non-standard scaling law 

or a breakdown of a standard scaling law, on the other hand, can often lead to the finding of a 
new universality class in physical systems. Recently, we found that a statistical quantity related 
to fluctuations follows a non-standard scaling law with respect to system size in a synchronized 
state of globally coupled non-identical phase oscillators [Nishikawa et al., Chaos 22, 013133 (2012)]. 
However, it is still unclear how widely this non-standard scaling law is observed. In the present 

y\ . paper, we discuss the conditions required for the unusual scaling law in globally coupled oscillator 

systems, and we validate the conditions by numerical simulations of several different models. 
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I. INTRODUCTION 

One of the essential studies of physics is to reveal universal phenomena in nature |l| . Synchronization is a universal 
concept in nonlinear science that is widely used to explain coherent behavior in a system of interacting elements [2| 
such as a reaction-diffusion system of chemical oscillators [3[ , an array of electrical units [J] , spontaneous luminescence 
of fireflies !5| , and a network of rhythmic biological components [6( • A synchronization transition in a system of N 
coupled oscillators takes place when the coupling strength surpasses a critical value at which an order parameter 
changes from zero to a non-zero value. Universality of synchronization phenomena can be partially confirmed by the 
fact that the scaling property of the order parameter in the vicinity of the critical synchronization transition point 
does not depend on the details of the system elements 0, l7H20l| . In the thermodynamic limit N — »■ oo, the study of 
the universal scaling law with respect to the coupling strength has been a topic of primary interest [3, Q ■ On the 
other hand, many studies have also focused on the universal scaling law with respect to finite system size N. In fact, 
it is known that the variance of the order parameter generally decays with N in the order of 1/N independently of the 
details of the system elements [9t4l2| . Such a scaling law is regarded as standard because it applies to a wide range of 
models. In the present paper, we focus on a non-standard scaling law that may provide an opportunity to determine 
a new universality class. 

To illustrate a difference between standard and non-standard scaling laws, we consider the globally coupled phase 
oscillator model UM, described as follows: 



K 



N 



-J2H6 k -6 j ){oTJ = l,...,N, 



N 



(1) 



fc=i 



where 6j represents the phase of the jth oscillator, loj represents the natural frequency of the jth oscillator, K > 
represents the coupling strength, and h denotes the coupling function. In the Kuramoto model where h(x) = sin(x) 
[lj,[l5|, the oscillators are synchronized when the coupling strength K is sufficiently large. Synchronization transitions 
have been extensively studied with regard to their analogy to the second-order phase transition [lj, LL5[ . A measure 
of synchrony is given by the order parameter R(t) as follows [14| |: 
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A synchronization transition is characterized by a change in the order parameter R(t) from almost zero to a finite 
value with an increase in the coupling strength K. In the thermodynamic limit N — > oo, the stationary state of 
R(t) in the incoherent (desynchronized) regime bifurcates at the critical coupling strength K = K c , above which 
the oscillators are coherent (synchronized). When N is finite, the effect of an increase in N on fluctuations can be 
different between the incoherent and coherent regimes [l|, [12] ■ To illustrate this, we introduce the diffusion coefficient 
of the time integral of R(t), characterizing long-term fluctuations of R(t) in systems of finite size [12j. We consider 
the variance of the time integral J Q R(s)ds, defined as follows: 



a 2 (t) = N 
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where (-) t represents the long-time average, and [-] s denotes the sample average over different realizations of uij that 
are chosen from a certain distribution. For the variance in Eq. ([3]), the following diffusion law holds (l2j ] : 



D = lim a 2 (t)/2t, 

t— VQO 



(4) 



which represents the deviation of L R(s)ds from its mean value {R)tt on the sample average. The standard scaling 
law of D with respect to system size N is described as follows [12| : 



D^ 0(1). 



(5) 



Namely, the value of D is almost constant independently of the system size N. We numerically confirmed that this 
scaling law is observed in the Ising model [l| , the probabilistic cellular automaton model [J, |2l], [22[ , and the random 
neural network model [23|. In the previous study [12 |. we also confirmed that the standard scaling law holds in the 
incoherent regime of the globally coupled phase oscillator model (jTJ . In our discussion in Ref . [12| , we suggested that 



the standard scaling law can be observed in a wide range of systems. We also found that, in the coherent regime of 
system ([1]), the asymptotic form of D for large N is represented as 

D ~ <3(l/iV Q ), (6) 

where a > is a certain constant. This non-standard scaling law suggests the existence of a new universality class. 
However, it is still obscure how broadly this scaling law is observed in coupled oscillator systems. The aim of this 
study is to elucidate the validity range of the above non-standard scaling law of D. 



II. NECESSARY CONDITION FOR NON-STANDARD SCALING LAW 

Let us consider the mean frequencies Uj — limt->oo(0j(t) — 9j(0))/t of the phase variables in the stationary state. 
We define by r(w) the distribution function of the mean frequencies of all the oscillators. A necessary condition for 
the scaling law ^ in system (JlJ is given as follows (12j |: 

lim r(w) = with N -> oo, (7) 

where $7 denotes the common frequency of the entrained oscillators. A typical example is illustrated in Figs. HJa) and 
[TJb), where we observe that Eq. (J7J holds in a coherent state whereas it does not hold in an incoherent state [2J. In 
Fig. [lja), oscillators with u ~ are dominant and oscillators with fast oscillations very few. In Fig. [ljb), we can 
observe the largest peak corresponding to entrained oscillators and two other peaks corresponding to non-entrained 
oscillators with fast oscillations. It has been previously shown that the value of D is mostly determined by the value 
of T(Q) around uj = Q 12]. The value of D usually decreases with system size N, if Eq. ([7]) holds. 

To understand the difference between the standard and non-standard scaling properties, we consider the following 
correlation function of the order parameter: 



C(t) = N 



Consequently, we can relate C(t) to D as follows (12ll25|: 



lim i I (R(t + t )-(R) t )(R(t )-{R) t )dt 
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Therefore, the form of the correlation function C(t) determines the value of D. In the incoherent regime of system 
(p}, the form of C(t) is almost exponential as shown in Fig. [TJc) [3, LL2[. In contrast, in the coherent regime, the 
non-entrained oscillators generate fluctuations in the correlation function C(t) around C(t) = as shown in Fig. HJd) 
pi LL2J- Therefore, a characteristic time period with C{t) < is observed. 



III. GENERALITY OF THE NON-STANDARD SCALING LAW 

Let us discuss what type of coupled oscillator models satisfy the condition ([7]). It has been reported that this 
condition is not satisfied in a coupled oscillator model with local couplin g I18II . one with very non- uniform coupling 
strength [24j, and one with a large variance in the natural frequencies [27|. To avoid these cases, it appears to 
be necessary that a model involves global coupling J2|, 0, |lj, [l6|, [30] , nearly uniform coupling strength [31| , and a 
sufficiently small variance in the natural frequencies [29j. 

We examine the existing models that are known to satisfy the condition |J7J in order to confirm whether the scaling 
law ([6]) holds or not. The results of our numerical simulations, presented in the following part, suggest that the 
condition (UJ) is not only a necessary but also a sufficient condition for the non-standard scaling law. 

First, we investigate the phase oscillator model with time-delayed coupling [28|], described as follows: 

K N 
fe=i 



where r represents the time delay. Equation J7J holds for any time delay r (28|. Second, we consider the phase 
oscillator model with inertia |26[ given by 

K N 
mBj + 6j =Uj + jjJ2 sin (^ " e i)> ( n ) 

fe=i 

where m corresponds to the strength of the inertia. The inertia introduces an adaptive effect in the system [26]. 
Numerical simulations have confirmed that Eq. ([7]) holds even for large values of m, i.e., for strong adaptive effects 
26] . Third, we explore the effect of amplitudes by using the globally coupled Stuart-Landau equations [27] , described 
as follows: 

K N 
Zj = (1 - \zj\ 2 + iu^Zj + — ^(z fe - Zj), (12) 

fc=i 

where Zj € C represents the complex state variable of the jth oscillator. The variance of the natural frequencies ujj is 
assumed to be not so large; otherwise, collective chaos or oscillation death occurs and violates Eq. ([7]) [27] . 

We show that all of these models satisfy the scaling law ([6]). The value of D is estimated by fitting a 2 (t) with a line 
of slope 2Dt except for the transient period, as shown in Fig. [5] For all the systems considered in the present paper, 
we clearly see that D is scaled asD~ 0(l/iV a ), as shown in Figs.[3Ja)-(c). In addition, the correlation function has 
a characteristic time period in which C(t) < as shown in Figs. [3Jd)-(f). Therefore, it is suggested that Eq. ([7]) is a 
sufficient condition for D ~ 0(1/ N a ) in globally coupled oscillator systems with large N. 

IV. DISCUSSION 

We have investigated the condition under which the non-standard scaling law ^ for the statistical quantity D 
holds. We have considered the relation between the non-standard scaling law and the condition in Eq. (J7]), and we 
have performed numerical simulations of several models satisfying the condition. The numerical results suggest that 
the condition is not only a necessary but also a sufficient condition for the non-standard scaling law. Our non-standard 
scaling law © does not sensitively depend on details of coupling schemes [l2| as well as on other factors such as time 
delay, inertia, and amplitude variables as demonstrated in this study. It is noteworthy that such a non-standard 
scaling law that can be observed in a large class of globally coupled oscillator models has not been reported thus 
far. In fact, even for conventional statistical quantities such as the variance and the correlation time of the order 
parameter, scaling laws with respect to coupling strength K and those with respect to system size N have not been 
fully explored except for those of the phase oscillator model dJ) 0, l9Wl2l. Ha . fl7H20 |. Thus, our result is useful for 
understanding fluctuations in general systems of globally coupled oscillators. 
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FIG. 1. Numerical results of T(tli) (upper panels) and the time evolution of the correlation function C(t) (lower panels) for 
the Kuramoto model (Eq. JTJ) with h(x) — sin(x)) with N = 64000. The natural frequencies uij are randomly drawn from a 
Gaussian distribution with mean zero and standard deviation one. (a) and (c) An incoherent state with K — 0.8. (b) and (d) 
A coherent state with K = 1.63 




FIG. 2. The time evolution of the variance a 2 (t) of the integrated order parameter in the coherent regime of system (|12[) with 
amplitude variables, where N = 8000 and K — 0.31. The value of a (t) increases linearly with slope 2Dt after a transient 
period. The natural frequencies ujj are randomly generated from a Gaussian distribution G{ui) with mean and standard 
deviation 0.2. The red line shows an average over 300 different initial conditions. 
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FIG. 3. Behavior of statistical quantities in the coherent regime. The upper panels show the scaling properties of D with respect 
to the system size TV ranging from 8000 to 64000. The lower panels show the correlation function C(t) for N — 64000. Each 
plot is an average over 100 different initial conditions. The diffusion coefficient D is scaled as D ~ 0(1/N a ) with a ~ 0.940 
for (a), a ~ 1.012 for (b), and a ~ 0.927 for (c). (a) and (d) System (|10p with time delayed coupling, where time delay r = 1 
and K — 1.78. We used the order parameter given by Eq. ([2}. (b) and (e) System (|ll|l with inertia term, where m — 0.1 
and K — 1.93. We used the order parameter given by Eq. @. (c) and (f) System (|12[l with amplitude variables, where 



K — 0.325. We used the order parameter defined as R(t) = ^ . 



/N. The natural frequencies loj are randomly generated 



from a Gaussian distribution G(u) with mean zero and a fixed standard deviation. The standard deviation of the natural 
frequencies uij is set at 1 for systems (|10[) and Ijllj l, and at 0.2 for system (|12|l . In the case where the natural frequencies are 
given deterministically i.e. j/(N + 1) = f^ 3 G(Cj)dCj, the result does not change. 



